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Abstract 

By treating idealized zeolite frameworks as periodic mechanical trusses, we show that the number 
of flexible folding mechanisms in zeolite frameworks is strongly peaked at the minimum density 
end of their flexibility window. 25 of the 197 known zeolite frameworks exhibit an extensive 
flexibility, where the number of unique mechanisms increases linearly with the volume when long 
wavelength mechanisms are included. Extensively flexible frameworks therefore have a maximum 
in configurational entropy, as large crystals, at their lowest density. Most real zeolites do not 
exhibit extensive flexibility, suggesting that surface and edge mechanisms are important, likely 
during the nucleation and growth stage. The prevalence of flexibility in real zeolites suggests that, 
in addition to low framework energy, it is an important criterion when searching large databases 
of hypothetical zeolites for potentially useful realizable structures. 
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Zeolites are crystalline aluminosilicate materials that are microporous to small molecules 
Because of their importance in the petrochemical and fine-chemical industries as cata- 
lysts and molecular sieves, there are ongoing efforts to synthesize new zeolite materials 2)]. 
At present there are 197 recognized zeolite framework types jjj, and that number has been 
doubling every 14 years. Recently, computer-generated databases of hypothetical zeolite 
frameworks have appeared with many millions of topologies with low calculated framework 
energies (46]. Given this abundance of energetically plausible structures, it has been puz- 
zling why only a small subset of these structures actually occur in nature. 

Real zeolites have complex structures and compositions. However, to explore their flexi- 
bility they can be idealized as tetrahedrally-coordinated SiO"2 frameworks comprising perfect 
Si(0*4)i/2 tetrahedra. Oxygen atoms are shared between tetrahedra and form the linkages. 
In real zeolites, tetrahedra are indeed close to perfectly regular. The intertetrahedral Si-O-Si 
bond angle is weakly constrained near 145°, but with a wide distribution of allowable angles. 
Dove and colleagues 5, Q showed that if tetrahedral silicates are modeled with force-free 
spherical joints at the oxygens, zero-frequency rigid-unit modes (RUMS) arise. Sartbaeva 
et al. [9 1 showed that many RUM modes correspond to finite folding mechanisms, allowing 
the stress-free framework to explore a range of densities termed the flexibility window. 



Remarkably, almost all of the known aluminosilicate zeolites exhibit a flexibilit y w indow 
[lo| . and the structures of real zeolites adopt the low density end of the window jillioj]. This 



state minimizes the framework electrostatic energy, largely by repulsion between oxygen 
anions [11 ]_. However, fewer than 10% of the low energy hypothetical zeolites in one large 
database 5| are flexible when represented as pure silicates. 

Tetrahedral silicate frameworks are locally isostatic, since the six of degrees of freedom 
per tetrahedron are exactly matched by the six constraints on the oxygens (three constraints 
per oxygen, which are shared by two tetrahedra). The flexibility in zeolites arises mainly 
from their high crystallographic symmetry, which imposes redundant constraints; in terms 
of rigidity theory, zeolites have non-generic structures [si, 12]. Periodicity adds six unit 
cell degrees of freedom [k]]. Three of these correspond to rigid-body translations. The 
remaining three, however, do not correspond to rigid-body rotations, as such motions are 
incommensurate with the lattice 13j. These correspond to three internal flexing degrees 



of freedom. Consequently, all periodic frameworks in three dimensions are guaranteed at 
least three infinitesimal folding mechanisms [13I [l^] , regardless of additional crystallographic 



2 



symmetry. These correspond to the three acoustic modes 15J . Infinite aperiodic frameworks 



are not guaranteed these mechanisms. An account of flexibility of periodic frameworks with 



internal symmetry has been given recently by Malestein and Theran [16| and Ross et al. (17|. 

To explore the configurational space of the flexibility window we employ techniques from 
rigidity theory by regarding a zeolite framework as a periodic truss. The 0-0 edges of 
tetrahedra are treated as bars of ideal length L = 0.263 nm, and the six-coordinated oxygen 



atoms behave as force-free spherical joints 
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18]. Si atoms are 'carried' inside the rigid 



tetrahedra and are ignored. The compatibility matrix, C, describes the 3N rigid constraints 
that maintain the bar lengths with N oxygen atoms per unit cell. The vector of infinitesimal 
displacements of the oxygen atoms, d, is related to the vector of 0-0 edge extensions, e, via 

Cd = e. (1) 

The condition e = describes a rigid framework. C is a 3iV x 3N matrix whose elements are 
Xij = (xj — Xi)/L, i/ij, Zij etc, with appropriate adjustments for periodic continuity. If we 
admit unit cell degrees of freedom, d acquires six additional components representing small 



changes in the unit cell dimensions, 5a, 5b, 5c, 5a, 5/3 and £7 [131 ] . C is now a 3iV x (3iV + 6) 
matrix embodying 3N + 6 degrees of freedom and 3N constraints. With compatible initial 
coordinates that satisfy e = 0, Cd = can be interpreted as an eigenvector problem with 
eigenvalues equal to zero. C is rank deficient by 6, yielding a null space containing at least 
six null solutions d (mechanisms), three of which will be the trivial translations 13(. The 
null eigenvectors are the modes of the framework that do not strain the 0-0 linkages. As 
the framework folds, the e = condition holds to first order between each step, but C must 
be continually updated when finite mechanisms explore their range. 

When there is no applied applied tension field (zero pressure) , the equation of motion is 

md = -kC T Cd. (2) 

m is the oxygen mass (we ignore Si atoms). We ignore Si-O-Si bending forces, so k is 
the spring constant of the 0-0 'bars'. Solutions of the type dj(i) = d; exp(ia^t) yield the 
eigenvalue equation 

C T Cd, = ^0,^ = 1^, (3) 

which is valid for small displacements. Removing the trivial translations, C T C is a (3N + 
3) x (3iV+3) square matrix. This underconstrained eigenvalue equation can be solved using 



singular value decomposition (SVD) [18J, and the null modes, oji = 0, define the infinitesimal 
mechanisms d; consistent with the flexibility window, e = 0. 

Thermodynamically, heat provides the energy that allows the system to explore its allowed 
configurations. In the limit of high temperature T, the vibrational entropy of a system of 
uncorrelated oscillators is 

3N+3 

(4) 



The product of eigenvalues equals the determinant, so the vibrational entropy is then 



S- 



ilog( 7 37V+3 det[C T C])|, (5) 
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where 7 = (h/k-QT) 2 k/m. This represents an entropic density, since (jSJ) applies to each 
configuration C within the flexibility window. Although this expression includes contri- 
butions from vibrational modes, the oui — null modes will dominate the entropic density 
since the logarithm diverges. To avoid divergences in (jSJ), we added 10~ 5 to the null modes 
(ui), rendering them 'floppy'. For large-amplitude mechanisms, we do not expect (jSJ) to be 
accurate. 

To explore the dependence on the size of the periodic unit, the unit cell can be expanded 
along the three axes to generate a supercell containing N c unit cells. The number of modes 
increases linearly with N c . Depending on the framework, the folding mechanisms will grow 
as iV™/ 3 (0 < n < 3). When n = 3 the entropic contribution from null modes is extensive, 
growing linearly with volume. The computational effort generally increases with the cube 
of the volume explored, i.e. ~ N%. To speed up computation, we applied Bloch's theorem 
d[r + R) = d(r) exp(2,7riK ■ R) to a single unit cell to explore the long wavelength null 
mode states within the Brillouin zone, and not just those at K = |.19l]. R is the position 
vector of a neighboring cell. The wavevector K has components that are multiples of 
(l/(pa),l/(qb),l/(rc)), where p, q and r are the cell repeats along the a, b and c axes 
(thus N c = pqr). Because of loss of periodicity, Bloch's theorem can not be used to follow a 
mechanism to finite amplitude. We use it to identify mechanisms, but follow them as modes 
at the T point (K = 0) in the p, q, r-expanded cell. 

Our computational tool, Zeolite Nullspace Explorer (ZeNuSpEx), follows the mechanisms 
as the structure evolves through the flexibility window. To initialize C, geometrically- relaxed 
initial coordinates were obtained using the GASP program [2^]. GASP binds the oxygen 
atoms to the vertices of ideal tetrahedra by springs of natural length, and then minimizes 



TABLE I: The twenty-five extensively flexible zeolite frameworks. The number of mechanisms at a 
general reciprocal lattice point is given. The '*' indicates that the framework could not be relaxed 
with perfect tetrahedra. 



ACO 
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EMT - 
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LTA 
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OFF - 
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SSF - 
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AST - 


12 


ERI - 
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LTL 
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PAU - 


24* 


TSC - 


25 
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FAU - 


16 


MER - 


2 


RHO 


4 


UFI - 


2 


DFO 


4 


KFI - 


6 


MOZ 


2 


SAS - 


2 


UOS - 


2 


EAB 


1 


LEV - 


3 


MTN - 


2* 


SOD - 


2 


uoz 
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the spring energy using the conjugate gradient method. C was adjusted continually by 
applying second-order corrections as the structure evolves. Occasionally, C is re gula rized 
by applying GASP. The singular value decomposition is the rate-determining step 2l|. We 
applied ZeNuSpEx to the 197 known zeolite frameworks and found that 25 of them exhibit 
extensive flexibility (Table H]). 

As an example, Fig. [T^i shows the SOD framework in its maximum symmetry (im3m), 
minimum density state, 16.58 T atoms/nm 3 . Figs. [T]d and QJ show the framework folded 
along two K = mechanisms to their densest configurations, 27.98 and 22.54 T atoms/nm 3 , 
respectively. The minimum density is determined by the onset of stretching of 0-0 edges; 
the maximum density is determined by oxygen overlaps, assuming an oxygen radius of 
L/2 — 0.135 nm. The body-centered unit cell of SOD supports 10 finite internal mechanisms 
at K = 0. This agrees with Hammonds et al. if we add their V and H modes for the primitive 
cell 15|. The minimum density state (Fig. [T^l) represents a confluence for all mechanisms, 



and so the density of states for null modes is highest at this end of the flexibility window. 
Few modes reach the densest states. 

As we increase the supercell size (K ^ 0) so that it contains N c primitive cells, the 
number of mechanisms increases linearly with system volume. A mechanism for a 4 x 3 x 1 
supercell of SOD is depicted in Fig. [T]i. All vibration modes and mechanisms can be 
excited simultaneously at small amplitude at the maximum unit cell volume condition. 
Once a mechanism is excited at large amplitude, lowering the symmetry and increasing 
the framework density, mechanisms that are incompatible with the lowered symmetry are 
suppressed. 
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Monte Carlo simulations of random mechanism excitations in SOD allow the system to 
explore the flexibility window. With increasing iV c , the system spends significantly more time 
near the minimum density state since fewer random moves will correspond to cooperative 
folding modes; the system spends most of its time where the density of mechanisms is 
highest. 

The contribution to the configurational entropy can be assessed by plotting log(z/j) /N$i 
against the mode index, sorted in order of increasing frequency (Fig. [2]). In Fig. [2^ we 
show the K = mode distribution for the minimum density and a folded configuration in 
PI symmetry. The folded structure, having PI symmetry, retains just 3 null modes. The 
shaded region represents the loss in configurational entropy when the structure folds. The 
equivalent plot for a 2x2x2 unit cell repeat (Fig. [2b) shows that the total number of modes 
increases to 52. There is a decrease in the shaded area, but with further increases in system 
size the mode distribution remains similar to that shown in Fig. [2b. In the asymptotic limit 
of large N c , the count of mechanisms grows extensively, two mechanisms per unit cell. 

Extensive flexibility is found to be associated with structures that have high symmetry. 
Conversely, PI guarantees just three infinitesimal modes, regardless of system size. Frame- 
works with higher topological symmetry always exhibit the highest density of mechanisms 
at the low density end of the flexibility window. A topologically high-symmetry structure, 
once folded to PI, is guaranteed only three mechanisms, of which at least one must be finite 
since this folded state is accessible from the high symmetry state. 

As a rule, the minimum density configuration is also the maximum symmetry state. 
The highest symmetry generally has the largest degeneracy of constraints, and thus the 
excess of degrees of freedom over constraints is maximum. Every mechanism supported by 
the structure has a path that passes through the minimum density state, which therefore 
represents the confluence of all mechanisms. Mechanisms can be linearly superimposed, 
with infinitesimal amplitude, at this lowest density state. Not all mechanisms are capable 
of reaching the high density end of the flexibility window, and so there is a gradient of null 
states that peaks at minimum density, and usually decreases monotonically as framework 
density increases. Fig. [3] shows how the entropic density term, S(p), depends on framework 
density, p, near the low-density end of the flexibility window. The curve is the average 
over 7000 different paths through the window. Superimposed is the average number of 
mechanisms available P(p) at each density. Both curves show a smooth decrease away from 
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the minimum density state. Earlier, Duxbury et al. 22j reasoned that the negative of the 
number of floppy modes in Bethe lattices behaves as a free energy. As we show, in zeolites 
extensive floppy modes reduce free energy by increasing the configurational entropy. 

Almost all of the 197 known zeolite frameworks have a flexibility window, whereas 12% 
of them are extensively flexible (see Tabled]). It appears that the existence of a flexibility 
window alone is important, and that the extensivity of flexible modes, despite an apparent 
thermodynamic advantage of increased configurational entropy in bulk crystals, is not essen- 
tial. As pointed out by Hammonds et al. 12j local flexible deformations can accommodate 
cations, or structure-directing molecules. Similarly, the surface and edge mechanisms of the 
flexibility window, which grow non-extensively as N^ 3 and N^' 3 respectively, may play a 
crucial role at the nucleation and growth stage of zeolite synthesis where N c is small. In 
our database of hypothetical zeolites, only about 10% of the low framework energy frame- 
works exhibit a flexibility window. It is clear that a flexibility window is the hallmark of a 
viable zeolite framework, and provides an important selection criterion in addition to low 
framework energy when evaluating databases for realizable structures. 
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FIG. 1: (a) SOD cage in its minimum density (and maximum symmetry Im3m) flexible state, (b) 
A higher density folded state that preserves /43m symmetry, (c) A folded state in PI symmetry 
at its maximum density, (d) A system of 4x3x 1 unit cells of SOD, viewed down [001], showing 
a higher-order mechanism in PI symmetry. The spheres are the centers of oxygen atoms. 
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FIG. 2: Plots of log(fi) per Si04 tetrahedron for the SOD framework, sorted by increasing fre- 
quency, for the minimum density configuration (solid lines) and a randomly-selected folded state 
(dashed lines), (a) Single unit cell (K = modes). There are 75 modes, 10 of which are null modes 
at maximum symmetry, I3m, reducing to 3 in the folded PI symmetry, (b) 2x2x2 supercell of 
the SOD framework. 579 modes, 52 null modes, reducing to 3 in the folded state. 
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FIG. 3: Plot of entropic density S(p) as a function of p. S(p) is sharply peaked at the minimum 
density p = 16.61 T atoms/nm 3 , where the framework symmetry is maximum. All mechanisms 
pass through this state. Superimposed on the plot is the frequency of states P(p)) arising when 
7000 random-walk paths are followed throughout the flexibility window. The flexibility window 
has one mode that reaches the highest density of 30.29 T atoms/nm 3 . The random paths never 
reached this state, but frequently encountered the minimum density state. 
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